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Abstract. In this article we study the pointwise decay properties of solutions to the 
wave equation on a class of stationary asymptotically flat backgrounds in three space 
dimensions. Under the assumption that uniform energy bounds and a weak form of 
local energy decay hold forward in time we establish a t~ 3 local uniform decay rate 
for linear waves. This work was motivated by open problems concerning decay rates 
for linear waves on Schwarzschild and Kerr backgrounds, where such a decay rate has 
been conjectured by Price 1391 . Our results apply to both of these cases. 



1. Introduction 

In this article we consider the question of pointwise decay for solutions to the wave 
equation on time independent, asymptotically flat backgrounds. Our interest in this 
problem comes from general relativity, more precisely the wave equation on Schwarzschild 
and Kerr backgrounds. There the conjectured decay, see Price |39] , is i -3 , and numerous 
partial results have been obtained. Our main result applies to both the Schwarzschild and 
Kerr space-times, and proves a stronger form of Price's conjecture, where sharp bounds 
are provided in a full forward light cone. However, the result does not depend on the 
specific geometry of either Schwarzschild and Kerr space-times; instead, it applies to a 
large class of stationary asymptotically flat space-times. 

We use (x,t = xq) for the coordinates in R 3+1 . We use Latin indices i,j = 1,2,3 for 
spatial summation and Greek indices a,/3 = 0, 1,2,3 for space-time summation. In R 3 
we also use polar coordinates x — rui with w e S 2 . By (r) wc denote a smooth radial 
function which agrees with r for large r and satisfies (r) > 2. Wc consider a partition of 
R 3 into the dyadic sets A k = {(r) « 2 fc } for k > 0. 

All operators we consider in this paper have time independent coefficients. To describe 
their spatial size and regularity we use the standard symbol classes S(r k ) = S : 

a G S(r k ) \dia(x)\ < c» fe ^, j > 

By S ra d{r k ) we denote spherically symmetric functions in S(r k ). We also introduce a 
slightly stronger class Z 1 S(r* ; ) with some dyadic summability added, 

a e ^(r fc ) ^2 TO ^- fe )||^a(x)|Uoo (Am) < Cj , j > 
m 

On occasion we allow for logarithmic modifications, and use the class of functions ^(log r) : 
aGS(logr)^=H a |<log(r), |^o(x)| < c 3 {r)~\ j>l 
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We consider the pointwise decay problem for the forward wave equation 

(1.1) (n g + V)u = 0, u(0)=«o, d t u(0) = Ul 

associated to d'Alembertian O g corresponding to to a Lorenzian metric g, and to a 
potential V. We work with potentials V of the form 

(1.2) V = V lr +V sr , V lr € S rad (r~ 3 ), V„€l 1 S(r~ a ) 
For the metric g we consider two cases: 

Case A: g is a smooth Lorenzian metric in R x R 3 , with the following properties: 

(i) g is time independent, i.e. dt is a Killing vector field. 

(ii) The level sets t = const are space-like. 

(iii) g is asymptotically flat in the following sense: 

g = m + g sr + gir, 

where m stands for the Minkowski metric, gi r is a long range spherically symmetric 
component, with S ra &{r ) coefficients, of the form 

gir = gir,tt(r)dt 2 + gi r ,tr(r)dtdr + gi^ rr (r)dr 2 + gi r ^(r)r 2 dui 2 

and 

(J ST IS cL short range component, with I Sir 1 ) coefficients, of the form 
9sr = g s r,ttdt + 2g 

We remark that these assumptions guarantee that dt is time-like near spatial infinity, 
but not necessarily in a compact set. This leads us to the second case we consider: 

Case B: g is a smooth Lorenzian metric g in an exterior domain K x K 3 \ 5(0, Ro) 
which satisfies (i),(ii),(iii) above in its domain, and in addition 

(iv) the lateral boundary R x 8B(0, Ro) is outgoing space-like. 

This latter condition insures that the corresponding wave equation is well-posed for- 
ward in time. This assumption is satisfied in the case of the Schwarzschild and Kerr 
metrics in suitable advanced time coordinates, such as the ones used in [30], [49]. There 
the parameter R is chosen so that < i?o < 2M in the case of the Schwarzschild metric, 
respectively r~~ < R < r + in the case of Kerr (see [TO] for the definition of r^), so that 
the exterior of the R a ball contains a neighbourhood of the event horizon. 

We consider the inhomogeneous forward Cauchy problem for the wave equation with 
initial data at time t = 0, 

(1.3) (a g + V)u = f, u(0) = uo, d t u(0)= Ul 

We denote the Cauchy data at time t by u[t] = (u(t),d t u(t)). In the sequel H k = H k (R 3 ) 
is the usual Sobolev space, and by H k ' 1 we denote the functions u with X7 x ii G H k . For 
the above evolution we begin with the following definition: 
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Definition 1.1. We say that the evolution (|1.3[) is forward bounded if the following 
estimates hold: 

(1.4) |H*]|| 

It suffices to have this property when / = 0. Then the / term can be added in by 
the Duhamcl formula. The uniform forward bounds above are easy to establish provided 
that the Killing vector field dt is everywhere time-like. Otherwise, there is no general 
result, but various cases have been studied on a case by case basis. 

In the case of the Schwarzschild space-time, following earlier work of Laba-Soffer |27j . 
Bluc-Soffer [3], [1], Twainy |52j and others, full uniform energy bounds were obtained by 
Dafcrmos-Rodnianski |16| in the exterior region. For a shorter alternate proof, as well 
as a discussion concerning the extension of the energy bounds to a neighbourhood of the 
event horizon we refer the reader to the article [3D] of the author and collaborators. 

As far as the Kerr space-time is concerned, partial results were first obtained by 
Finster-Smoller |18j . and Finster-Kamran-Smoller-Yau [17] . Full uniform energy bounds 
were proved in the the case of the Kerr space-time with small angular momentum in 
Dafermos-Rodnianski [2]; their result applies in effect to a larger class of small pertur- 
bations of the Schwarzschild space-time. An alternate proof of the Kerr result, based on 
local energy decay estimates, was given shortly afterward in Tataru-Tohaneanu [35] , 

A stronger property of the wave flow is local energy decay. We introduce the local 
energy norm LE by 

(!-5) \\ u \\le = sup||(r)~5 U || L2(RxAfc) 

k 

its H 1 counterpart, 

\\u\\ LE i = \\Vu\\ LE + \\(r)-\\\ LE 

as well as the dual norm 

(i-6) II/IUb* =ElK r )'/IMHx^) 

k 

These definitions are specific to the three dimensional problems. Some appropriate mod- 
ifications are needed in other dimensions, see for instance |34j . 

We also define similar norms for higher regularity, 

"Ir- 
respectively 

\\f\\LE^ = E W^fhE' 
\a\<k 

The dual spatial norms C£ and C£* are the spatial counterparts of the space-time LE 
and LE* norms, 

(1.7) \\v\\c £ = S up\\(r)-h,\\ LHAk) \\f\\a>. =J2\\(r)?f\\mA k ) 

k k 
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Their higher regularity versions are defined as 

rn m 

(1.8) IMUen. = £sup||<r>-*ajt;|| iW ||/|k^ = ^ElK r ) l ^/H^) 

J=0 fe j=0 fe 

In Case A above this leads to the following 

Definition 1.2. a) We say that the evolution (|1.3[) /ias £/ie local energy decay property 
if the following estimate holds: 

(1-9) IMUf^ <c fc (||u[0]|| ff i. fcxfffc + ||/|| LB ..0, fc>o 

mRx K 3 . 

The first local energy decay estimates for the wave equation were proved in the work 
of Morawetz [35] , [36] , [37] . There is by now an extensive literature devoted to this topic 
and its applications; without being exhaustive we mention Strauss [43], Keel-Smith- 
Sogge [22], [23], Burq-Planchon-Stalker-Tahvildar-Zadeh [6], Metcalfe-Sogge [31] . [32] . 

The sharp form of the estimates as well as the notations above are from Metcalfc- 
Tataru [34] ; this paper also contains a proof of the local energy decay estimates for small 
(time dependent) long range perturbations of the Minkowski space-time, and further 
references. See also the related paper Metcalfe- Tataru [33], as well as Alinhac's article pQ. 

There is a related family of local energy estimates for the Schrddingcr equation. As a 
starting point we refer the reader to [48] and references therein. 

Work in progress of Metcalfe- Tataru suggests that in the nontrapping stationary case 
the local energy decay (jl.9p is closely tied to a non-resonant spectral behavior at frequency 
(see [29] for our similar results for Schrddingcr evolutions). 

In Case B an estimate such as (|1.9[) is unlikely to hold due to the existence of trapped 
rays, i.e. null geodesies confined to a compact spatial region. However a weaker form 
of the local energy decay may still hold if the trapped null geodesic are hyperbolic. 
This is the case for both the Schwarzschild metric and for the Kerr metric with angular 
momentum \a\ < M. To state such bounds we introduce a weaker version of the local 
energy decay norm 

\Hlei = 11(1 - x)^Ale + IKr^ttlUB 
for some spatial cutoff function x which is smooth and compactly supported. Heuristi- 
cally, x is chosen so that it equals 1 in a neighbourhood of the trapped set. We define as 
well a dual type norm 

\\f\\LE^ = \\xf\\L^m+Y,W( r ^fhHA k ) 
fe 

As before we define the higher norms LE^ k respectively LE^; k . 

Definition 1.3. We say that the evolution (|1.3[) has the weak local energy decay property 
if the following estimate holds: 

(1.10) \HleIi* ^ c k (\\u[0}\\ Al , kxHk + \\f\\ LK ,*), k > 

in either IxR 3 or in the exterior domain case. 
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Two examples where weak local energy decay is known to hold are the Schwarzschild 
space-time, and the Kerr space-time with small angular momentum \a\ <C M. 

The estimates in the Schwarzschild case are the end result of a series of papers by Laba- 
Soffer [27], Blue-Soffer [3], @],[33], Blue-Sterbenz [5] Dafermos-Rodnianski [T5], [T3] and 
Marzuola- Metcalfe- Tataru-Tohancanu [3U] . In particular we note the contribution of [TB] 
to the understanding of the red shift effect near the event horizon, as well as that of 
[38], |30] concerning improved bounds near the photon sphere (i.e. the trapped set). The 
sharpest form of the estimates is obtained in |30| ; in particular the bounds there fit into 
the framework of Definition 11.31 

In the case of the Kerr space-time with small angular momentum \a\ <C M the local 
energy estimates were first proved in Tataru-Tohaneanu |49j , in a form which is compati- 
ble with Definition 11.31 A more detailed presentation as well as stronger bounds near the 
trapped set are contained in Tohaneanu's Ph.D Thesis [50] as well as in the follow-up pa- 
per [51] . For related subsequent work we also refer the reader to Dafermos-Rodnianski [12] 
and Andersson-Blue [2J. 

The high frequency analysis of the dynamics near the hyperbolic trapped orbits has 
a life of its own, but not much to do with the present article. For more information 
we refer the reader to work of de Verdiere-Parisse [11], Christianson [9], Nonnenmacher- 
Zworski [52] and references therein. 

We remark that the local energy decay property is stronger than the mere boundedness 
of the evolution. This is not as clear in the case of the weak local energy decay; however, 
one case which is well understood is when the Killing vector field dt is time-like within 
the support of the cutoff function \ m t ne definition of LE^. 

In this article we consider the pointwisc decay of solutions to the wave equation (|1.3[) . 
If g is the Minkowski space-time then, by Huygens principle, if the initial data has 
compact support then u will vanish in a compact set after some time. If instead we add 
a stationary potential which has exponential decay at infinity and no eigenvalues or zero 
resonances, then one obtains exponential decay of solutions inside any compact set. For 
polynomially decaying stationary perturbations of the Minkowski space-time the question 
of the precise decay rates is still largely open. In principle, one expects that these decay 
rates will be determined by three factors: (i) the spectral behavior near time frequency 
(e.g. eigenvalues, resonances), (ii) the trapping properties of the metric and (iii) the 
behavior of the metric at spatial infinity. 

This work was inspired by an idea in Tataru [48] and Metcalfe- Tataru [34], which is 
roughly that the local energy estimates contain all the important local and low frequency 
information concerning the flow, and that only leaves the analysis near spatial infinity 
to be understood. Precisely, in [35] and [33] it is proved that local energy decay implies 
Strichartz estimates in the asymptotically flat setting, first for the Schrodinger equation 
and then for the wave equation. The same idea was exploited in |30] and [51j to prove 
Strichartz estimates for the wave equation on the Schwarzschild and then on the Kerr 
space-time. 
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The result we obtain in this article shows that the above philosophy applies as well 
to the pointwise local decay. Precisely, we prove that if weak local energy decay esti- 
mates hold, then the pointwise decay rate of the solutions to the wave equation is only 
determined by the behavior of the perturbation at spatial infinity. 

To define the initial data spaces we use several vector fields: the generators of spatial 
translations T = {81,82,83}, the generators of rotations ft = {xidj — Xjdi] and the 
spatial scaling S r = rd r . Then we define the weighted Sobolev spaces Z m ' k with norms 

(1.11) ||/||z».»= sup \\{r) n ™S k r f\\ C£ * 

i-\-j+k<rn 

Now we are ready to state the main result of this article: 

Theorem 1.4. Let m be a large enough integer. Let g be a metric which satisfies the 
conditions (i), (ii), (Hi) mix R 3 , or (i), (ii), (Hi), (iv) in R x R 3 \ 5(0, Ro). Let V be 
a potential as in (jl.2[) . Assume that the evolution (|1.3|) is forward bounded and satisfies 
the weak local energy decay estimates (jl.101) . Then in normalized coordinates the solution 
u to the homogeneous problem (jl.ll) satisfies the bounds 

(1.12) KM)I < (f + N) ^_ N)a (1Kb— + NIL-) 

respectively 

(1.13) IM*,*)I S (t + N) |,_ N)3 (1Kb— + 1Kb-) 

The normalized coordinates referred to in the theorem are introduced in the next 
section via a spherically symmetric change of coordinates. These are needed in order to 
insure that the outgoing light cones are close to Minkowski cones. 

The required decay of the initial data at infinity is roughly a summable r -3 for the 
position uq, respectively a summable r for the velocity u\. These are sharp. By 
contrast, we made no effort to optimize the value of m. 

In view of the local energy decay estimates in Marzuola-Metcalfe-Tataru-Tohaneanu 
[50] . respectively Tataru-Tohaneanu [IH], the hypothesis of the theorem is valid in the 
case of the Schwarzschild space time, and also the Kerr space-time with small angular 
momentum, with V = 0. Thus we obtain: 

Corollary 1.5 (Price's Law). The decay estimates (|1.12[) . respectively (|1.13j) hold for the 
Schwarzschild space-time, as well as for the Kerr space-time with small angular momen- 
tum, with respect to coordinates which coincide with the Regge- Wheeler coordinates near 
spatial infinity and with the Eddington-Finkelstein coordinates near the event horizon. 

This problem has had a long history. Dispersive L 1 — > L°° estimates providing t^ 1 
decay of 3 + 1 dimensional waves in the Minkowski setting have been known for a long 
time. There is an extensive literature concerning dispersive estimates for □ + V; this is 
by now a well understood problem, for which we refer the reader to the survey paper of 
Schlag [£Q. 

The need for weighted decay inside the cone arose in John's proof [5D] of the Strauss 
conjecture in 3 + 1 dimensions. Decay bounds for □ + V , similar to those in the theorem, 



LOCAL DECAY OF WAVES ON ASYMPTOTICALLY FLAT STATIONARY SPACE-TIMES 7 

were obtained by Strauss-Tsutaya [H] and Szpak [JB], [35] i see a l so Szpak-Bizoh-Chmaj- 
Rostworowski [37] . 

A more robust way of proving pointwise estimates via L 2 bounds and Sobolcv inequal- 
ities was introduced in the work of Klaincrman, who developed the so-called vector field 
method, see for instance |24| . This idea turned out to have a myriad of applications, and 
played a key role in the Christodoulou-Klaincrman |10j proof of the asymptotic stability 
of the Minkowski space time for the vacuum Einstein equations. 

In the context of the Schwarzschild space-time, Price conjectured the t~ 3 decay rate 
for linear waves. More precise heuristic computations were carried out later by Ching- 
Leung-Sucn- Young [7J , [5] • Following work of Wald [53] , the first rigorous proof of the 
boundcdness of the solutions to the wave equation was given in Kay- Wald [21] . 

Uniform pointwise t decay estimates were obtained by Bluc-Stcrbenz |38j and also 
Dafcrmos-Rodnianski [TS] ; the bounds in the latter paper are stronger in that they extend 
uniformly up to the event horizon. Very recently, a local t~ s decay result was obtained by 
Luk [28| . These results are obtained using multiplier techniques, related to Klainerman's 
vector field method; in particular the conformal multiplier plays a key role. 

Another venue which was explored was to use the spherical symmetry in order to 
produce an expansion into spherical modes, and to study the corresponding ode. This 
was pursued by Kronthaler [26], [25], who in the latter article was able to establish the 
sharp Price's Law in the spherically symmetric case. A related analysis was carried out 
later by Donningcr-Schlag-Soffcr [40] for all the spherical modes; they establish a t~ 2 
local decay for each spherical mode. 

Switching to Kerr, the first decay results there were obtained by Finster-Kamran- 
Smoller-Yau [IT] . Later Dafcrmos-Rodnianski [T^] were able to extend their Schwarzschild 
results to Kerr, obtaining almost at -1 decay. 

Finally, we mention the related work of Dafermos-Rodnianski [15], where Price's Law 
is established in a nonlinear but spherically symmetric context. 

We conclude the introduction with a brief overview of the paper. In the next section 
we discuss the normalization of the coordinate system, which is done in order to insure 
that the outgoing null cones are close to Minkowski cones. In the following section we use 
the uniform forward energy bounds in order to define resolvent operators R T associated 
to the evolution (jl.lj) in the lower half-space {3r < 0}, and to establish straightforward 
L 2 and H m bounds for the resolvent. These bounds degenerate as 3r approaches 0, so 
they provide no information about the resolvent on the real axis. 

In Section |4j we transfer the local energy decay estimates to the resolvent R T . These 
bounds have the key advantage that they are uniform as r approaches the real axis. 

Using a version of the vector field method, in Proposition 15.101 we commute the re- 
solvent with the rotations f2 and scaling S to obtain a larger class of resolvent bounds. 
These bounds allow us in Proposition 15 . 1 II to prove that we can extend the resolvent R T 
to the real axis continuously in a weaker topology as well as show that it inherits all the 
estimates previously proved in the lower half space. We also establish the key outgoing 
radiation condition (|5.34[) for R T for real r. 
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The limit Rq of the resolvent R T as r — !• is considered in Propositions 16. 13116.141 this 
allows us to define Rq, and establish some key properties for it. The precise description 
of Rq is used in Proposition 16.151 to produce an expansion up to second order for the 
resolvent R T near r — 0. 

In Section [7] we switch from L 2 to L°° bounds in the resolvent estimates for real r, 
in a way which is reminiscent of the Klaincrman-Sobolev embeddings. Exploiting the 
outgoing radiation condition (|5.34[) . we also obtain improved decay estimates at spatial 
infinity for (d r — ir) m R T . 

In the last section we conclude the argument. Assuming for simplicity that u = 0, the 
time Fourier transform of the solution u to the forward wave equation (|1.1[) is represented 
as u = R T u\. Thus in order to invert the Fourier transform and obtain pointwise estimates 
for u we need to control the size and regularity with respect to r of R T u\ . This is relatively 
easy for large r because we allow the initial data to have high Sobolev regularity in 
Theorem 11.41 The delicate issue is the resolvent behavior at r = 0. In order to obtain a 
tr 3 decay for u at infinity, we use the second order expansion for R T near r = given 
by Proposition 16. 1 51 In particular we are able to identify the leading contribution of the 
long range part of the wave operator and of the potential. 



2. The coordinate normalization 



Here we carry out a preliminary step which brings the metric to a canonical form near 
spatial infinity and removes certain logarithmic type corrections in the resolvent bounds 
later on. This is achieved in two steps, a spherically symmetric coordinate change fol- 
lowed by a conjugation with respect to a scalar weight. This transformation is roughly 
equivalent to using Rcggc- Wheeler coordinates in Schwarzschild/Kcrr near spatial infin- 
ity. 

Our final goal is to reduce the problem to the case when the operator D g + V in (jl.l[) 
is replaced by an L 2 selfadjoint operator P which has the form 

(2.14) P = -d 2 + A x + Pi r + P sr 
where the long range spherically symmetric part has the form 

(2.15) P lr = gt{r)^ + VJ r (r), V lr € S rad { r - 3 ) 
and the selfadjoint short range part P sr has the form 

(2.16) P a r = d a fi?dp+V ar , g a J efSir- 1 ), Vsre^Sir- 3 ) 
with = 0. 

2.1. Normalized coordinates. We first remark that the condition (hi) translates into 
a similar representation for the dual coefficients g a ° , 

We can simplify these coefficients via a radial change of coordinates as follows: 

Lemma 2.6. There exists a spherically symmetric change of coordinates so that the long 
range component g"f of g a ^ has the additional properties 

(2.17) gf r = -gj r r , gg = 
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for large r. 

We will refer to the coordinates given by this lemma as normalized coordinates. These 
are the coordinates used in Theorem 11.41 In Schwarzschild or Kerr space-times these 
coordinates are precisely the Regge- Wheeler coordinates for large r. In a bounded spatial 
region, any choice of coordinates which is nonsingular on the event horizon will suffice. 



Proof. To achieve the relation g\l = we keep r unchanged and reset t by setting 

t:=t + a(r) 

where a is determined by 

rt 

o'(r)=- : 9lr 



1 + 

Then we have a' £ £rad(?" _1 ) and a € S ra d[hir), therefore the property (iii) is left 
unchanged. We remark that this can only be done provided that 1 + gf^ > 0. This is 
clearly true for large r, and can be arranged for small r by modifying the short range 
part P sr . 

Once g\l = 0, we insure that g\^ = with a radial change of coordinates 

dr * = [ ] dr 



.i + C 

This gives r* — r € 5 rat ;(lnr) and dr* /dr € 1 + -Srac*^ -1 )- Thus while both transforma- 
tions above leave condition (iii) unchanged, both the functions a(r) and r*{r) — r may 
have a logarithmic component as r —¥ oo. The proof is concluded. 

□ 



2.2. Left and right multiplication by scalar functions. Consider a metric g as in 
(12.171) . The corresponding d'Alembertian has the form 

\9 

where it is easy to see that g E 1 + S r ad(r ) + ^S^r -1 ). This is turned into an L 2 
selfadjoint operator via the conjugation 

U g -> g?n g g-i 

To further insure that the coefficient gf r of d 2 is set to —1 we multiply the above conju- 
gated operator by (— g lt )~ l E 1 + ^S^ 1 ); in order to preserve the L 2 selfadjointness, 
half of it will go on the left and half will go on the right. 

Thus, our final replacement for D g + V is the operator 

P = 9k-9 tt )-H° 9 + V)(-g tt )- h 9-* 

Commuting, P can be written in divergence form, and it will have the required form 
given by QZJfy . (g/EjJ) and (|2TT5]| . 

Also it is easy to see that such a transformation leaves unchanged the three properties 
of the forward evolution (|1.3[) : boundedness, local energy decay and pointwise decay 
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3. The resolvent and energy estimates 

Consider the solution u to the forward homogeneous problem 

(3.18) Pu = 0, u(0) = 0, u t (0) =g € L 2 , 

extended by to negative times. 

By the uniform energy bounds its time Fourier transform u(r) is a H 1 valued distri- 
bution in R, and admits a holomorphic extension into the lower half plane. We define 
the resolvent operator R T by 

R T g = u(t), 3r < 
The uniform energy bounds (|1.9p for P translate into an L 2 bound for R T , 

(3-19) WRrgWfr. + t\\Rt9\\l* < TcC-rllfflU'. »r < 0, 

\^ST\ 

as well as similar bounds for higher Sobolev norms, 

(3.20) \\Rrg\\H^+r\\Rrg\\^<T^\\g\\H^ %T<0. 

In order to write an equation for the resolvent we express the operator P in the form 

P = -d 2 + P l d t + P 2 

where, in view of (|2.14p . (|2.15|) and ()2 . 16[) . P 1 is a short range spatial skew-adjoint 
operator and P^ a second order spatial operator of the form 

P 2 = A + P lr + P 2 r 

with both long and short range components. 

Since u is extended by to negative times, it solves the distributional equation 

Pu = g5 t =o 

Taking a Fourier transform in this relation, it follows that u(t) solves the equation 

(3.21) P t u(t) = g, P t = t 2 - irP 1 + P 2 
Hence for the resolvent we obtain the equation 

(3.22) v = R T g P T v = g 

By Duhamcl's formula we also obtain a representation for the Fourier transform of the 
solution v to the full forward Cauchy problem 

(3.23) Pv = f, v(0) = h, v t (0)=g 

Lemma 3.7. Let v be the solution to (|3.23l) with f G L°°L 2 , h 6 H 1 and g 6 L 2 . Then 
for 9r < we have 

(3.24) v(r) = RtU{t) +rh + 2P 1 h + g) 



This allows us to show that P T and R T are inverse operators if 3r < 0: 

Lemma 3.8. Let 3r < 0. Then P T : H 2 — > L 2 is a bounded one to one operator with 
dense range, and R T is its inverse. 
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Proof. Given u € H 2 , apply the previous lemma to the function u = u lt>o to sec that 
R t P t uq = uq. The range is dense since it contains H 1 . Indeed, by (|3.22[) and (|3.20[) , for 
/ G H 1 wc have R T f G H 2 and P T R T f = f. □ 

4. The resolvent and local energy decay 

The goal of this section is to repeat the arguments of the previous section, but using 
the local energy decay property instead of the uniform boundedncss. For this we consider 
the solution u to the inhomogcneous forward problem 

(4.25) Pu = ge LE*> k+i , k>0 

with g supported away from t = — oo. From the local energy decay estimate (|1.10[) we 
obtain the weaker bound 

(4-26) MIls 1 -* < c fc||/||ii?*.fc+3 

with a loss of three derivatives. We will use this bound to derive a similar resolvent bound. 
For each r we use the C£ k norms defined in (| 1 . T[) . (|1.8[) to introduce the r dependent 
norm C£ T by 

(4.27) IMI^j = \\{\t\ + {r)- l )u\\ C£k + \\Vu\\ C£k + \\(\r\ + (r^J-Vu^ 
Then we have 

Proposition 4.9. Let 3r < 0, k > 0, g G C£*' k+A and v = R T g. If the weak local energy 
decay bound (|4.26[) holds then we also have 

(4.28) \\v\\c£i<ck\\g\\c£*^ 



While this estimate is for now stated for 3t < 0, it is nevertheless uniform as r 
approaches the real axis. Later we will use this fact to extend this estimate to all real r. 
If — > 1 then this bound is superseded by (|3.20[) . Thus only the case — 1 < 3t < is 
of interest. 

Proof. If we consider / supported in a fixed dyadic region Ai and measure u in another 
region A m , from (|4.26p we obtain the uniform dyadic estimates 

J2 2-^\\\7^u\\ LHMxAm) +2'^\\\7l t u\\ L 2 {RxAm) < J2 ^WKJhHR^) 

i<k 2<fc+3 

Since the equation (|4. 25[) is solved forward in time, it is straightforward to add an ex- 
ponentially decreasing weight in the above estimate. Precisely, from the above estimate 
one obtains a weighted version 

^2-^|| e - £ *V^ 1 U |U 2(KxAm) +2-^|| e - £i V^ U || L2(R><Am) < ]T 2^||e- et V^/|j i2(Rx ^ ) 

i<k i</c+3 

which holds uniformly with respect to e > 0. If wc further restrict e to < e < 1 then 
we can uniformly move the exponential inside the differentiation to obtain 

^2-^||V;+ 1 e^|U 2(KxAm) +2-^||V^e-^|| L2(K><Am) <^ 2^ || V^e- e */|| i2(ExAi) 

i<k i<k-\-3 
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Since we only have Hilbcrt space norms in the above estimate, we can take a Fourier 
transform and use Plancherel's theorem to obtain 

2-f ^(1 + |r|) fc -' (||Vi +1 M(r - ie)\\ L2{RxAm) + (2~ m + \r\)\\%u(r - ie)\\ L2(WLxAm) ) 

i<k 

< 2Hl + \r\) k - l \\V l J(r- l e)\\ L2{RxAl) 

i<k+3 

where the L 2 norms are now taken with respect to t 6 1 and x £ K 3 . 

We recall that by (|3.24[) we have u(t — ie) = R T f(r — ie). Observing that the above 
bound holds for all / in a dense subset of the Hilbert space defined by the norm on the 
right, it follows that in effect its pointwise version with respect to r must also be valid, 
namely 

2"? £(1 + |r|) fc - (||V*>|| L2(Am) + (2— + \VtT\)\\Vl t v\\LHA m) ) 

i<k 

< 2Hl + \r\) k - l \\Vl t g\\ LHAl) 

i<k+3 

Here we have rcdcnotcd r — ie by r, which now satisfies 3r < 0, and replaced /(r) by g 
and u(t) by v = R T g. We remark that in view of (|3.20[) we can replace the factor |3?r| 
by |t| in the second term on the left. 

Taking the suppremum over all to > 0, and then summing over all / > 0, we arrive at 
(4.29) E^ + M^OI^II^ + IK^^ + H^M^ E ( 1 + \^) k+3 ~ i \\9\\ce^ 

This concludes the transition from local energy decay to resolvent bounds. Wc remark 
that at the formal level one can go directly from (|4.26[) to (|4.29|) for real r via Plancherel's 
theorem. However, in order to rigorously justify this analysis we arc using the interme- 
diate step of deriving uniform bounds for Qt < 0. 

It remains to prove that (|4.29|) implies (|4.28[) . This depends on the size of r. If |r| w 1 
there is nothing to do. In the two remaining cases this transition is carried out via simple 
elliptic arguments. 

(a) The case |r| <C 1. Then we need to estimate the second derivative of v for large r. 
Wc write the equation for v in the form 

P 2 v = -t 2 v + itP x v 

This allows us to estimate 

IK(r)- 1 + \r\)- 1 P 2 v\\ C£k < T\\v\\ C£k + HV^H^* + \\g\\ce*.k 

where all the right hand side terms are already controlled by the right hand side in (|4.26|) 
via (|4.29l) . Since P 2 is elliptic for large r, the transition from P^v bounds to V 2 i> bounds 
is done in a standard elliptic fashion within each dyadic region A m . 

(b) The case |r| ^ 1. The left hand side of (|4.29[) already contains the C£ k norm of 
v, we only need to be able to discard the r factors on the right. For this we observe that 
the equation for v, namely 

(t 2 + irP 1 +P 2 )v = g 
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is elliptic for low spatial frequencies |£| <C t. This suggests wc split g into a low frequency 
and a high frequency part, 

.9 = S< T (D x )g + S> T (D x )g := g tow + g high 

For the high frequency part we can estimate 

i</c+3 

and use (|4.29j) directly. 

For the low frequency part, on the other hand, we reiterate the equation, setting 



where v\ solves 



Vlow = t 2 giow + V\ 



P T vi = ir 1 P 1 giow + t 2 P 2 giow = giow,i 



Thus we have traded powers of r for derivatives, which is favourable since gi ow is low 
frequency. Reiterating several times we arrive at 

fc-i 

Vlow = E r 2 j Q]9low + Wl ow 
3=0 

where Qj are partial differential operators of order j with smooth bounded coefficients 
and wi ow solves 

PrWiow = T~ k Qkgiow + T~ k ~ 1 Qk+igiow ■= hi ow 
The first term in viow is estimated directly in terms of gi ow while for wi ow we remark that 

E (l + H) fc+3 - 1 ||/»w|| jC £M<||fl|| jC£ ..* + . 
i<k+3 

and use (|4.29j) . The proof of the proposition is complete. 

□ 

5. Vector field bounds and the resolvent limit on the real axis 

Our first goal in this section is to prove further estimates for the operator P T by 
commuting it with several vector fields: 

i) The rotations VL = {xidj — Xjdi}. Commuting them with P T we obtain 

\_Pt ; ^] Q sr 

where Q sr stands for a short range operator of the form 

Q sr = T(h M d t + d l h M ) + dih ij dj + h, h Q \ h ij e ^(r- 1 ), h G fS^' 3 ) 

ii) The scaling S = —t8 t + rd r . Then the commutator is 

[P T , S] = 2P T + Qir + Qsr 

where the long range component is radial and has the form 

Q lr = k"A u + k, k G S rad (r- 3 ) 
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The next proposition applies for spatial functions g, in which case S = S r g. However, 
it is also interesting to allow g to depend on the parameter r. 

Proposition 5.10. Let 3r < and g € C£ * , possibly depending on t, so that 

(5.30) \\T i n j S k g\\cs* < 1, i + Aj + 16k < N 
Then 

(5.31) \\T*WS k R T g\\c £T <l, i + Aj + 16k<N-4 



Proof. Set v = R T g. To illustrate the method we first consider the simplest cases. In the 
case of £1 we have 

P T Qv = Qg + Q sr v 

A direct computation shows that 

(5-32) \\QsA\ cs ** < \\v\\ C£ , 

Hence by (|4.28[) we obtain 

ll^llre* ^ i 1^.9 1 1 + IMLe-.fc+s 

Consider now the operator fl 2 . We have 

P T fl 2 v = il 2 f + Q sr nv + Q sr v 
so we conclude again via (|5.32[) and (|4.28p that 

W&vWcs* £ \\n 2 g\\ C£ ,, k+ , + ||n 5 || iC£ .,H-. + \\g\\ce?w 

Next consider the scaling field S. We have 

P T Sv = Sg + 2g + Q sr v + Qi r v 
For Q sr we use (|5.32|) . while for Qi r we use 

(5.33) IIQjrvll^..* < Mce* + l|n a «ll££* 
to obtain 

II^II^J < \\Sg\\c£*- k +^ + \\^ 2 g\\c£'-"+ s + H^IUf.fc+is + ||fif||££*.*+W 

The proof of the general case is by induction. The details are left for the reader. □ 

Now we have enough information in order to define the pointwise limit of the resolvent 
as t approaches the real line. The next result deals with the limit on 1\ {0}. The 
detailed analysis near frequency zero is left for the next section. 

Proposition 5.11. i) The operators R T extend continuously from the lower half-space to 
the real line M \ {0} in the H^ omp — > Lf oc topology, and the bound (|4.28|) holds uniformly 
for all real t. 

ii) Let t € M. \ {0} and g G ££*' 4 . Then the function v = R T g satisfies the outgoing 
radiation condition 

(5.34) lim 2-i\\(d r + iT)v\\ LHAj) =0 
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Hi) Conversely, suppose that r G R \ 0, and v € L£\ which satisfies the outgoing 
radiation condition (|5.34j) . If P T v = g € C£*' then v = R T g. 

This result allows us to transfer the previous bounds in this section from the lower 
half-plane to the real axis. 

Corollary 5.12. The results in Propositions ^. 9\5.10\ and Proposition \6. 15\ apply for all 
nonzero real t. In addition, for i,j,k as in Provosition \ 5.1(K the functions T 1 ^ S k R T g 
satisfy the outgoing radiation condition (|5.34|) . 

Proof, i) In view of the uniform bound (|4.28|) it suffices to establish the convergence on 
a dense subset. Let g g C£* so that (|5.30j) holds for some large N. Then by (|5.31[) we 
have 

\\SRrg\\ce T <l 

Hence we can write 

rd T R T g = SR T g + rd r R T 

In view of the L£ T bounds for R T g and SR T g we can bound the LHS locally in L 2 and 
in H 1 . Hence the map r — > R T g is locally Lipschitz away from r = from {3r < 0} into 
H} oc . The desired convergence follows. 

ii) By (|4.28j) it suffices to establish the radiation condition (|5.34|) for g in a dense 
subset of L£ *' . Here we consider functions g € ££* so that (|5.30[) holds for g for some 
large N, and so that in addition g decays one order faster at infinity, (r)g <G ££ *. We 
denote v(t) = R T g. 

For some large R we truncate w(r) to the exterior of a ball of size R and set 

w(t) = X>rv(t). 

This does not affect its behavior at spatial infinity and thus the radiation condition. The 
functions w(t) solve an equation of the form 

(5.35) (A + t 2 )w = Qi r w + Q sr w + [P T , x>r]v(t) + x>R9 ■= H T ) 

We can now use the local energy decay estimate for A + r 2 , namely 
(5-36) \\Vw\\c£ + \\(\t\ + (ry^wWcs < \\h\\ L]S . 

which can be either proved directly or derived from (|1.9|) for □ in the same way (|4.28l) 
was obtained from (jl.IOl) . 

In view of (|5.31|) and of the condition (r)g € C£ * we can estimate uniformly for r in 
a compact subset of {3r < 0} \ {0} 

\\{r)h{T)\\ C s* < 1. 

The radiation condition follows if we can show that the solutions w(t) to (|5.35j) satisfy 
(5-37) \\{r){d r +iT)w\\c£<\\{r)h\\ C£ * 

uniformly on compact sets. Indeed, let us estimate w in a dyadic region A m . We split 
h = X<m-2h + x>m.~2h. We have 

\\x>m-2h\\ce- <2- m \\{r)h\\ ce , 
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therefore by (|5.36|) we have 

\\Xm(dr + ir)(A + r 2 )- 1 x> m -2^||££ < 2- m \\(r)h\\ C£ * 

It remains to prove the same bound for the other component x<m-2h. For this use 
the kernel for the operator (A + t 2 )" 1 which is 

K T (y) = \v\- X e~ iT ^ 

This is justified for 5r < 0, where the symbol of A + t 2 is nonzero, and in effect bounded 
from below. A direct computation gives 

(d r +iT)(A + T 2 )- 1 h(x)^e- lTlxl J h(y)\x - y\- 1 e-^ x - y \a(x,y)dy 

where 

, . x-{x-y) . f x-(x-y) 

a{x,y) = -—- [J - it —- - 1 

\x\\x — y\ A \ \x\\x — y\ 

We are interested in the case where we have the localization « 2 m and 2\y\ < \x\. 
There we have \a\ < 2~ m (y). Taking absolute values inside the integral, we use again 
(|5.36|) with t = to obtain 

\\ Xm (d r + »t)(A + r 2 )- 1 x<m _ 2 / l || ££ < 2- m . 

The proof of (|5 . 37[) is concluded. 

iii) If v is compactly supported then for e > we have 

V = Rr- ie P T - ie V 

But the right hand side converges to R T g locally in L 2 , and the proof is concluded. 

Suppose now that v is supported away from 0. We first use the radiation condition 
for v to derive a similar radiation condition for its derivatives, 

(5.38) lim 2~i\\V k (d r + tT)v\\ L 2 [A]) = 0, fc<4 
We recall that P T has the form 

(5.39) P T = A + t 2 +P lr + P sr , 

therefore by commuting and estimating all terms with decaying coefficients in terms of v 
we obtain 

IKr)" 1 ^ + T 2 )(d r + it)v\\ c£ „ < \\v\\ C£i + \\g\\ c£ ^. 
Hence from the radiation condition (|5.34[) it follows that 

lim 2-i\\A(d r + iT)v\\ L 2 {Aj) =0. 

j-s-oo 

Combined with (|5.34j) and elliptic estimates for A we obtain a similar decay property for 
V(<9 r + it) and V 2 (9 r + it). Reiterating we arrive at (|5.38|) . 

For e > we define the functions v e = we _e ^ r ^ which satisfy v t = R T ^i e P T ^i e v e . By 
(|4.28p and part (i) of the proposition, the desired conclusion would follow if we are able 
to show that 



(5.40) Pr-tcVe -> 9 in C£ 



*,4 
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For this we use (|5.39[) to compute directly 

P T - le v e = ge-^ r) +2e{d r + iT+-)ve-^ r) + {el 1 S(r- 2 ) + e 2 l 1 S{r- 1 ) + el 1 S(r- 1 )V)ve- t( > r) 

r 

Since g G ££*' 4 , the first term converges to g in C£*' 4 . For the (d r + ir)v part of the 
second term, the same follows from the expanded radiation condition (|5.38|) . Finally for 
the remaining terms we obtain 0(e) decay in ££*' from the ££ T bound for v. Thus 
(|5.40p follows and the proof is concluded. □ 

6. The resolvent analysis near zero frequency 

In this section we first study the limit of the resolvent at frequency r = 0. Then 
we consider data with more regularity and decay at spatial infinity, and we produce an 
expansion of the zero resolvent in terms of powers of r. Finally, we obtain a quadratic 
expansion of the resolvent near frequency for data with faster decay. We begin by 
establishing the existence of the limit of the resolvent at zero. 

Proposition 6.13 (Resolvent extension to r = 0). a) Let g,(r)g G ££*' 4 . Then the 
limit 

R g = lim R- ie g 
exists in the Lf oc topology, and the following bound holds: 
(6-41) \\(r}Rog\\c £ * < \\(r)g\\£ £ *A + *, k>0 

b) The operator Rq admits a unique continuous extension to ££ *' 4 which satisfies 
(6-42) \\Ro9h £ * < \\9\\c £ *A +k , fc>0 

respectively 

(6.43) lim 2^-^\\V m R Q g\\ LHA]) = 0, j = 0,1,2 

j->oo 

c) Conversely, let v G C£q so that Pqv = g € ££*' 4 and 

(6.44) lim 2-*||0 r u|| L 2 (A . ) = 0. 
Then v = R^g. 

Proof, (a), the case k = 0. Denote v e = R-i e g. By the local energy decay estimate 
(|4.28[) we know that 

This already suffices in order to obtain the existence of a strong limit v e — > v in H} oc on 
a subsequence. To prove convergence it suffices to establish the uniqueness of the limit. 
For this we first study the regularity of the limit. A-priori we know that v G ££q, but 
this docs not suffice for uniqueness. 

The equation for v e has the form P_i e v e = g, or in expanded form 

(P 2 - eP 1 - e> e = g 
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with P 2 sclfadjoint and P 1 skew-adjoint. Since v e £ H 2 , we cam multiply by v e and 
integrate by parts to obtain 

5? J gv f dx = J —P 2 v € ■ v € + e 2 v 2 dx 

The principal part of — P2 is positive definite outside a compact set K, so we obtain 

liv^lll. < 

For the first term on the right we use the fact that (r)g £ C£* and (r) _1 w £ £ C£. For 
the second term on the right we use the C£i e bound for v e . In the last term we use the 
Hardy inequality to absorb the far away part into the left hand side, and the C£i t norm 
for the near part. We obtain an uniform bound 

\\Vve\\h<\\(r)gh £ *,* 

so any limit v of a subsequence must belong to H 1 . Suppose we have two distinct limit 
points v\, V2- The difference v — v\ — V2 must satisfy v £ H , Pqv = 0. But this 
would contradict the local energy decay estimates for the function v lt>o which solves 
the homogeneous wave equation in [0, 00). 

It remains to boost the decay of v and its derivatives near infinity to (r)v £ C£q. For 
this we view Pq as a perturbation of —A near infinity. Precisely, for large R we consider 
the equation for w — x>R v > 

Pqw = X>n( r )g + [ p 0, X>R(r)}v := h 
Since the above commutator has compact support, we can estimate 

||<r>ft||£f < \\(r)g\\c£* + \\v\\c£„ < \\g\\ce--' 
We rewrite the above equation as 

(6.45) -Aw = h + X> R/2{r)(Pir+P S r)w 

We solve this last equation via a Picard iteration. It is easy to verify directly that 

IKrH-Ar^lU^IKr}/^. 

where (— A) -1 is the convolution operator with r _1 . 

On the other hand the last term on the right is perturbative for large R, 

\\(r)x>R/2(r)(Pir+P s r)w\\ C£ , < iT 1 1| (r)w\\ CSo 

Thus a unique solution w £ (r)~ 1 ££ C H 1 is found iteratively. On the other hand, by 
energy estimates as above the equation (|6.45j) also has a unique solution in H , therefore 
the (r) ~ l C£q coincides with the H 1 solution constructed from v. This proves the existence 
of the limit of u e as e — >• 0, as well as (|6.41[) for k = 0. 

(b), proof of (|6.42|) . The bound (|6.42|) is obtained by passing to the limit r = ie — > 
in the estimate (|4.28[) . Since (r) -1 ^*' 4 is dense in C£*' 4 , this provides us with an unique 
continuous extension which satisfies the bound (|6.42l) . 

(a), the case k > 0. Now (|6.41[) for k > 1 follows by combining the regularity of v 
provided by (|6.42[) for small r with the decay at infinity given by the k = case of (|6.4ip 



gv e dx 



\VvJ 2 dx 



K 



l dx 
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plus elliptic regularity for large r. The reason we need to argue in this roundabout way 
is that Pq is elliptic for large r, but not necessarily for small r. 

(b) , proof of (|6.43[) . By (|6.42[) . it suffices to establish (|6.43|) for g in a dense subset 
of C£* A+k . But this is provided by (|OT1) . 

(c) . This proof is identical to the proof of Proposition 15.111 (iii). First the higher 
regularity version of (|6.44l) are established, namely 

lim 2-i\\V k d r v\\ LHAj) = 0, fc<4. 

Then the relation v = P^g is obtained as the limit of v e = R^.i e P-i e v e for v e = ve~ tr . □ 

Our next result provides better bounds for Rq on functions with more regularity and 
decay. 

Proposition 6.14 (Higher regularity for R ). a) If g £ Z k+iA then the following repre- 
sentation holds for large r 

(6.46) R g = c(r)(r)- 1 + Vl 

where, with N arbitrarily large, we have the bounds 

(6-47) + ||<5'rc||;is(i) + l^illz^- 1 ^ !l.9llz fc + 4 > 1 

b) If we further have g E z k+4 ' 2 then Rof admits the large r representation 

(6.48) R g = c(r)" 1 + d(r) • V(r)- 1 + e(r)(r)- 2 + v 2 
where, with N arbitrarily large, we have the bounds 

(6.49) |c| + ||d||ioo + \\Srd\\iiso.) + ||e|| S (i) + W^Wz"." < hWz^ 

Here, by a slight abuse of notation, by the S(l) and Z 1 5(l) norm we understand the 
sum of an arbitrarily large but finite number of seminorms for the respective symbol 
classes. 

We note the key presence of the e(r)r -2 term in l|6.48[) . which is the highest order 
contribution due to the long range scalar potential. This term will later turn out to 
produce the first nonanalytic contribution in the expansion of R T with respect to t near 
t = 0. Thus this term ultimately determines the final decay rate in our main theorem. 

Proof, a) Consider (|6.46[) . For the local (small r) regularity of v we only need the bound 
(I4.28P . Thus it suffices to prove the same estimate for the solution w to the near infinity 
problem (|6.45[) . We begin with the simpler case h € (r)~ l C£* , in which case we seek to 
prove that w admits a representation as in (|6.46[) so that 

(6.50) ||c|U~ + IM^) +^||(r)-WV> 1 |U £ . < \\{r)h\\ C£ , 

J'<2 

The second right hand side term in (|6.45[) is perturbative with respect to this estimate 
with a small 0(R^ 1 ) bound, so it suffices to prove that this representation is valid for 
the simpler equation 

-Aw = h 
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Then w is computed using the fundamental solution for —A in R 3 , namely 

w = h * r _1 

We consider a dyadic decomposition of h with respect to the {.A m } partition of' 

h = y^h m 

m 

Correspondingly wc partition w into 



For the first part we directly estimate in a scale invariant elliptic fashion 
2- 2 ™\\w l ° w h> +2- m \\Vw l ° w \\ L , + UVVTIb < IlUb 
Hence is easy to see that the sum of the w l ™ contributions can be included in V\. 
Consider now the functions w^l gh , which we expand near infinity as 

(6.51) U)'™ 9 ' 1 = C m X>m+2\x\~ 1 + Z m, 

where 

c m = J h m (y)dy, z m = x>m+2 J h m (y) ( J~~^ ~ ]^|) dy ' 

By Cauchy-Schwarz we can estimate 

IMIi* £ \\{r)h\\ce' 
therefore the first term in (J6.51[) yields the first term in (|6.46l) . with 

C ( r ) = £ C '"X>rn+2 
m 

The second term in (|6.51[) . on the other hand, can be incorporated into v±. Indeed we 
can bound 



1 1 



\x-y\ \x\ 
which yields the off-diagonal decay 

\\(ry^z m \\ Laa{An) < 2" l ~", m>n + 2 

which suffices after summation with respect to n and m. This concludes the proof of 
(|6.50p . The more general bound (|6.47[) follows from (|6.50|) by elliptic regularity once we 
observe that -A(c(r)r~ 1 ) e /^(r -2 ). 

b) Finally we consider the case when g E Z k+i ' 2 . Again it suffices to prove the same 
estimate for the solution w to the near infinity problem (|6.45[) . In a manner similar to 
the previous case, we will first construct a decomposition (|6.48[) of w so that the following 
estimate holds: 
(6.52) 

| C | + ||d||^ + ||^|| ;ls(1) + J R^||e|| s(1) +Ell« 2 ^V> 2 ||^ <\\h\\ {r) - 2C£ , +R i s d(r _ 4) 

j<2 

Once this is done, its higher regularity version (|6.49[) follows by elliptic theory. 

We remark that in the above estimate we have harmlessly added an S ra d{f~ A ) com- 
ponent to h. Its purpose is to render the second right hand side term in (|6.45[) fully 



LOCAL DECAY OF WAVES ON ASYMPTOTICALLY FLAT STATIONARY SPACE-TIMES 21 

perturbative. The R^i weights serve the same goal. To understand how this works we 
split w into w = c(r) + w\, and consider all possible terms: 

i) The output of both the short range operator P sr and the long range operator Pi r 
applied to w\ can be estimated as follows: 

\\X>R{P S r+Plr)™l\\ {r) -*C£* < R' 1 (II 4 L™ + II S r d\\ ,i s(1) + || e\\ + ]T || (r) 2 ^ V> 2 j| ) 

which provides an gain compared with (|6.52[) . 

ii) The output of the short range operator P sr applied to c(r)(r) _1 satisfies the weaker 
bound 

\\PsrC(r) - 1 \\p 8 (r-V < 

which gives decay when restricted to r > R, 

|| P^r)- 1 1| (r) - 2££ . <o R (l)\c\, 

where the decay of the constant as R — > oo comes from the fact that the I 1 summation 
is restricted to dyadic regions A m with 2 m > R. 

iii) The output of the fc' J (r)A aj part of Pi r applied to c(r)(r) -1 is zero. 

iv) Finally he output of the k(r) part of P; r applied to c(r)(r) _1 belongs to S(r~ 4 ), 
and we have an R~z gain compared to (|6.52[) due to the choice of constants in (|6.52l) . 

Thus we need to prove the bound (|6.52l) for an expansion (|6.49[) of the solution w to 
the equation 

—Aw = h, supp h C {r > R} 

We consider first the case when (r) 2 h £ ££* . We proceed as in the proof of (16.461) - 
(16.471) . The analysis of w l ™ rests unchanged. However, for w% 9h we need a second order 
expansion at infinity, namely 

(6.53) w^ 9h = c m x>m+2\x\~ 1 + dl nX>m +2d l r~ 1 + z m , 



where 



h m (y)dy, d l m = J y t h m (y)dy 



Then we estimate 



m|| L£* 



2 m \c m \ + \d) n \ + \\z m \\ ZN , a <\\(r) 2 h. 
Hence we can set 

c = ^c TO , d\r) = s ^d l m x> m +2(r), e(r) = r^c m x< m +2{r) 

m rn m 

One easily checks that c, d l and e satisfy bounds as in (|6.52[) . so the proof of (|6.52[) is 
concluded in this case. 

It remains to consider the case when h € S ra d{r~ A )- Then w is also smooth and 
spherically symmetric, and we can write 

d 2 (rw) = rh 
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Integrating this relation twice, first from infinity and the second time from zero, it follows 
that w has the form 

w 2 ^c(r)- 1 +e(r)(r)" 2 , e £ S ra d(l) 
Furthermore, since h is supported in r > R, we obtain 

|c| < R- 1 

Thus the proof of (|6.52[) for w is concluded. □ 

The next step is to use the Rq estimates in order to improve the analysis of R T g for r 
near 0, and in particular the bounds in Proposition 15.101 in the case when g has better 
decay at infinity. Precisely, as in the previous proposition, we will successively consider 
the case when g E z k+i ^ and when g G Z k+i ' 2 . We will use the notation a A b for a 
smooth 1-homogeneous function of a, b G R + which equals min{a, b} unless a ~ b. 

Proposition 6.15. a) Let g € z k+4,1 and \t\ < 1 with 3r < 0. Then the function 
v(t) = R T g can be represented as 

v( T ) =e- lT ^v(0)+ Wl (T) 

where the function w\ is given by 

Wi{t) = R T (x>\T\-i(r)g + rh) 

with h satisfying 

(6.54) \\r' n {d r + iT) m T l VL j S l h\\ C £> < 1, m + i + j + I < k 

b) Assume in addition that g € Z k+4,2 . Then v can be represented as 
v(t) = (v(0) + TV! + Te (r,T))e~ lTr + w 2 {t) 

where 

V!=R gi, gi£Z k '\ 
the radial function eo has the form 

eo(r, r) = r^e^r A Irl" 1 ) + r(e 2 (r A It^ 1 ) - e 2 {\T\- 1 )) 1 e 1 ,e 2 £ S(logr) 
and u>2 is given by 

w 2 (t) = i? r (x>|T|- 1 .9 + T 2 h) 

with h satisfying (|6.54|) . 

Proof, a) We note that the function e^ lT ^v(0) belongs to the range of R T since it satisfies 
the outgoing radiation condition (|5.34|) . We compute the equation for the function wi = 
v(t) -e- 2Tr v{0): 

P rWl = g + (ge-^) - g ) + 2r(d r + l)v(0)e-^) 

(6.55) r 

+ (Tl l S{r- 2 ) + r 2 / 1 ^^- 1 ) + r/^r-^V) v{Q)e~ lT ^ 

where the second term on the right comes from the Laplacian and the third is generated 
by the short range part of P. We observe that the long range part of P does not produce 
any contributions because the function e lT ^ is radial. 
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For v(Q) we use the representation (|6.46|) and the corresponding bounds (|6.47p . This 
allows us to estimate all terms on the right but the first one as in (|6.54|) . We note in 
particular that in order to estimate the contribution of the first term in (|6.46p to the 
second term above, it is necessary to take advantage of the improved ^5(1) bound on 
S r c in (1Q71) . 

Finally we consider the first term on the right in (|6.55[) . which we rewrite as 

9e- ir {r) - 9 = X<\r\-ig(e- tT{r) - 1) + X>\T\-ige- lT{r) - X>\t\-i9 

In the first term we can pull out a r(r) factor from the difference 1 — e~ %Tl<r \ and then 
(16.541) follows. For the second factor, (|6.54p is obtained directly. 

b) We start again with (|6.55|) . but now we can only place into h the terms which have 
a t 2 factor; we call such terms negligible errors. This includes the middle term in the 
last expression in (|6 . 55[) . Now v is as in f|6 .48[) - (|6 .49[) . therefore the remaining part of the 
same expression can be written in the form 

T{l x S{r- 2 ) + / 1 ^(r- 1 )V)w(0)e- ir < r > := Tgi e- lTi - r \ g x g Z k ' x 

and use the argument in part (a). Since g\ already has a factor of r, the error in part (a) 
will have a factor of t 2 . We further note that, compared to part (a), the error Tt/ie _lT ^ 
is better because it already contains the oscillatory factor e~ lT ' r ' ; we call such terms type 
1 negligible errors. 

It remains to consider the second term in (|6.55[) . This is where we need to take 
advantage of the more precise representation for v(0) in (|6.48l) - (|6.49[) . The contribution 
of cr _1 vanishes, while the contribution of v 2 in (|6.48p is type 1 negligible. 

Consider now the error term r(d r + (^di(r)dir^ 1 ) e~ lrl<r \ : arising from the second 
term in (|6.48p . We compensate for it with the explicit correction 

tv 2 = rdi (r)di (In r)e _iT < r > 
A direct computation shows that 

P T V2 - (d r + i) (diWftr- 1 ) e" iT < r > = (P T - (A + t 2 ))v 2 

+ drd^dir-^e-"^ + 9 2 d i a,(lnr)e~ 4T < r > 

It suffices to show that the right hand side can be expressed in the form <?2e~ tr ^ + rh 
with g 2 G l l S(r~ z ) and h as in (|6.54p . This follows by taking into account the effect of 
the short and long range terms in P T for the first term, and by using the bound on S r di 
in (|6.49p for the last two terms. 

Finally we consider the term r(d r + \:)e{r)e~ l ^ T € 5 f ~ 3 e~ ? ^ r . This is the term which 
arises due to the long range part of P T , and which is in effect responsible for the decay 
rate in our main theorem. Here there is no mechanism that can eliminate logarithmic 
type corrections for the resolvent. The redeeming feature of this case is that we can work 
fully with spherical symmetry, and that we can replace P T by A + t 2 modulo negligible 
perturbations. 

We begin by solving the equation which corresponds to r = 0, namely 

A l)3 = (cV + -)e(r) e S(r- 3 ) 
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and we obtain a solution of the form tj 3 = (r) ei(r) with e\ £ S ra d(logr). However, 
unlike all previous cases, here the function Tv^e~ lT ^ can no longer be a suitable correc- 
tion, because that would inherit the logr term even when r > r~ 1 . Instead we consider 
the modified function v 

v(r,r) = (r)~ 1 e 1 (r A |r| _1 ) 
and use v 3 = tv(t, r)e~ %T ^ as the correction. A direct computation shows that 

P T v 3 = T (d r + -)e(r)e~ lT{r) + 2r 2 {d r + -)ie~' lr{r) + Tg 3 (r) e - lT{r) + r 2 h 
r r 

Here g 3 <G l l S(r^) and h as in (|6. 54[) are negligible errors generated by the long range 
and short range part of P T . Thus only the second term provides non- negligible errors 
and needs further attention. Modulo terms which can be included in h, the expression 
(d r + \)v has the form 

(d r + -)v = x<M-i(r)(r)d r e 1 (r) + h 
r 

Even though this comes with a r 2 factor, it is not quite negligible because it only has 
decay (r)~ 1 d r ei G S(r~ 2 ) at infinity so it logarithmically fails to be in LE* uniformly 
with respect to r. Instead we need to reiterate again. 

We solve 

Ae 2 = (r)~ 1 d r ei 

and obtain a solution e 2 G S (log?'). Then we choose the last correction of the form 
v 4 = r 2 (r)- 1 (e 2 (r A |t|- X ) - e 2 (|r|- 1 )) e - ?;T(r) 

Here it makes no difference whether we add the factor e~ lTr or not, as U4 is supported in 
{ r l T l S 1} an d does not see any oscillations. Again, a direct computation shows that 

Prv 4 = r 2 X< | r |- 1 (r)(r)a- 1 e 1 e- lT < r > + r 2 h 

with h as in (|6.54p . a negligible error. Now the first term has a better r 2 factor, so it is 
also negligible. The proof of the proposition is concluded. 

□ 

7. From L 2 to pointwise bounds 

In this section we supplement the weighted L 2 resolvent bounds with corresponding 
pointwise bounds. Given a r dependent function g(r) which is defined in 3r < or in 
an open subset of it, we look for pointwise bounds for R T g and its derivatives. First we 
consider the case when g has only C£ * type decay at spatial infinity. 

Proposition 7.16. Let g £ ££* , depending on r , so that (|5.30l) holds, 
(i) Large t, \t\ > 1. Then the following pointwise bounds are valid: 

(7.56) \T i n j S k R T g\<(r}- 1 , i + 4j + 16k < N - 12 
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(ii) Small t, |t| < 1. Then the following pointwise bounds are valid, 

( mm{l, {\T\ir))- 1 } i = 

(7.57) \T i QPS k R T g\ < I (r)- 1 i = l i + Aj + 16k < N - 12 

[ (r)- 2 + |r|(r)- 1 i>2 

(Hi) In addition, if r is real then the outgoing radiation condition holds: 

(7.58) lim r{d r + iT)T l Vt J S k R T g = Q i + Aj + 16k < TV - 12 

\x | — too 

Proof, (i) Denote v = R T g and g ijk = T 1 WS k R r g, v ijk = T l WS k R T g. By Proposi- 
tion [5TTU1 g,ijk satisfy (|5.31[t . Using the Sobolev embcddings on the sphere 

II^IU°°(s 2 ) ^ II</>IIl2(s2) + II^ 2 0||l2(§2) 

we obtain for i + Aj + 16k < N - 12: 

(7.59) X^talU^SfCA™) < 1, i + Aj + 16k<N~& 

m 

respectively 

(7.60) ll^fe|U=L-(A m ) <2~^, i + 4j+16fc<iV-12 

Using these bounds we revisit the equation for Vij k which we can rewrite in the form 

(7.61) {d 2 + T 2 )(rv ljk ) = r^AuVijk + rQi r v<i,<j t < k + rQ ST .u<i,<j,<fc + rg<i,<j,< k 
In view of (|7.59[) and (|7.60l) we obtain 

(7.62) ^2T||(a r 2 + r 2 )(™ ufe )|| i 2 Ls3( ^ m) < 1, i + 4j + 16fc<iV-20 

m 

Using (|7.60|) for Vij k and d r Vij k we can localize Vij k so that the above bound is preserved, 

ll(<9 r 2 +T 2 ){xA m rv ljk )\\ L -2 L ^ (Am) < 
Now we use the following fundamental solution for d 2 — t 2 , 

K r {s) = t- 1 ^ 

This satisfies the bounds 

\K t {s)\<\t\-\ \d s K T {s)\<l 
Then the last bound leads to 

\ T \\\XA m rv ijk \\ L ^ {A:j) + \\d r (xA m rv ijk)\\L'*'{A j ) H 1 

which we rewrite in the form 

(7.63) M|i>«k| + |$.««fc|<<r>- 1 

This concludes the proof of (|7.56|) . 

(ii) The arguments above remain valid in dyadic regions A m with 2™ > | t - | 1 , and 
directly yield (|7.57|) for i = and i = 1; in the latter case we need the additional 
observation that T can be expressed as a bounded linear combination of d r and r _1 fi. 

Similarly, we can express T 2 in terms of d 2 , r~ 1 d r Q and r~ 2 Q 2 . Thus, in view of 
flTjTJl for i = 0, 1, for the i = 2 case of (f737)) it suffices to bound d?Q j S k v. For this 
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we return to (|7.61[) . move the term T 2 Vij k on the right and use (|7.57|) for i = 0,1. This 
allows us to obtain a favourable bound for all terms on the right except for the d 2 part 
of Q sr . We obtain 

\\d?v 0jk \\ Lao{Am) < \T\(r)- 1 +2- m \\d 2 r v ,< J ,< k \\ L ~ {Am) 

The second term on the right does not appear when j = and k = 0. Hence we can use 
induction with respect to j, k to obtain 

RWlU"(A m ) <\r\{r)^ 
and conclude the proof of (|7.57j) when i = 2. The case i > 2 is identical. 

It remains to consider the case when to is small, namely 2 m < |t| _1 . There the 
argument is slightly different. Again by the Sobolcv embeddings on the sphere, from 
(|5.31j) we obtain 

(7.64) 2- m \\v ijk \\ LlL oo (Am) + \\Vv ijk \\ L2rL ^ {Am) +2 m \\V 2 v ijk \\L^(A m ) < 2"? 

Then the bound (|7.57|) for i = 0, 1 follows by Sobolev embeddings with respect to the 
radial variable. Finally we obtain f|7. 5T[) for i > 2 directly from the equation (|7.6ip as 
above. 

iii) Suppose r is real. From the spherical Sobolcv embeddings and the L 2 radiation 
condition (|5.34[) applied to the functions Vij k we obtain 

(7.65) lim 2%\\(d r + iT)v ijk \\ L 2 L o«=0 

m— >oo T 

Writing d 2 + t 2 = (d r — ir)(d r + it) we interpret (|7.62[) as a first order ode for the 
functions (d r + ir)(rVij k ). The averaged bound (|7.65[) allows us to integrate this ode 
from infinity to obtain 

lim (d r + it) (rvijk ) = 

which implies (|7.58l) . 

□ 

Next we consider a more favourable case when g either has better decay at spatial 
infinity, or it has an oscillatory behavior which matches the one of the resolvent. For 
simplicity wc restrict ourselves to real r. If t is large then only the former alternative 
needs to be considered. In this case we will prove the following: 

Proposition 7.17. Let g € Z" hn with m n. Then for real r with \t\ > 1, the function 
v = R T g satisfies the following pointwise bounds: 

(7.66) \{Td T )\ve lT{r) )\ < \T\ l (r)- l -\ I < n 

Proof. We first observe that by Proposition 17. 16[ the bounds (|7.56j) hold with TV = to. 
Also, by Sobolcv embeddings, g satisfies the stronger pointwise bounds 

(7.67) \T i g\ < {r)- 2 -\ i<n-2 

For x in a compact set the exponential e* T ^ is harmless, and on the other hand we 
can write rd T as a smooth combination of S and tT. Hence the bound (|7.66|) follows 
from (|7.56[) . It remains to prove (|7.66l) for large r. 
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Applying (|7.56[) for i = 0, j = 0, k = we directly obtain (|7.66[) for I = 0. For I = 1 
we compute 

Td T (ve ZTr ) = Sve iTr + r(d r + ir)ve lTr 

The first term is estimated again by (|7.56|) . Hence we need to prove that for large r we 
have 

(7.68) \( dr + i T )v\<^ 

For ru we have in polar coordinates the equation 

(d 2 , + T 2 ){rv) = rg - r~ 1 A LJ v + rQi r v + rQ sr v 

Using ([736]) and ((737)) we can bound the RHS by |r|r" 2 . Then 

\(d r - ir)(d r + ir)(rv)\ < \r\r~ 2 

which is integrated from infinity using the radiation condition to obtain (|7.68[) . The same 
argument applied to the equation (|7.6f p for Vijk instead of v yields 

(7.69) | {d r + ir)v ljk \ < M { + 4j + lQk < m - 20 

For I = 2 we compute 

{Td T ) 2 {ve" r ) = S 2 ve lTr + r(d r + ir)Sve lTr + r 2 (d r + ir) 2 ve lTr 
The second term is handled via (|7.69[) so it remains to prove that 

(7.70) \(d r +iT) 2 v\<^ 
We have 

(d r - ir)(d r + ir) 2 (rv) = (d r + ir)(rg + r _1 A w i; + r~ 2 v + rQ sr v) 

= (g + r- 2 A^v + r~ 3 v + Q sr u) + r(d r + ir)f 

+ r~ 1 A L0 (d r + ir)v + r~ 2 (d r + ir)v + rQ sr (d r + ir)v 

Since (d r + ir)v also satisfies the radiation condition (|7.58p . it suffices to show that the 
RHS admits a r 2 r~ 3 bound. This follows from the pointwise bounds (|7.56p for u for the 
terms on the first line, and (|7.69p for the terms on the second line. 

A similar computation using (|7.61|) also shows that for large r we have 

(7.71) \(d r +iT) 2 v ijk \ < -j , i + 4j + 16fc<m-40 

The proof of (|7.66p is completed by induction with respect to /. □ 

Next we consider the case of small r. Here instead of the stronger decay assumption at 
spatial infinity we consider a weaker condition for g and establish two types of regularity 
bounds with respect to the spectral parameter r: 

Proposition 7.18. Let m be a large integer. Let g be a Z rn '° valued function of r in 
[—1, 1] which satisfies the additional bounds: 

(7.72) \\{r) n (d r + iT) n T l W S k g\\ C£ * < 1, n + i+j + k<m 
Then the function v = R T g has the following properties: 
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(i) r < |r| . Then we have 

(7.73) \{rd T ) l v\<l, Km 

(ii) r > \t\~ . Then we have 

(7.74) \{Td r ) l (ve lT ^)\ < (Irl(r))- 1 , I « m 

Proof. To prove the proposition we rely on the estimate (|7.57[) applied to the function v. 
For g on the other hand, by Sobolev embeddings we have the pointwise bound 

(7.75) \T i Q j S k (ge iT < r) )\ < (r)~ 2 ~\ i+j + k<m-2 

(i) . For I = the bound (|7.73|) follows directly from (|7.57|) . Consider now the case 
1 = 1. For small r we write 

S T v = — Sv + S r v 

which is again bounded by (|7.57|) . Even better, (|7.57[) shows that for the functions v^k 
we have 

(7.76) \S T Vijk\ + \S r Vijk\<l, i + 4j + 16k<m-2Q. 
For I — 2 we write 

S^Vijk = S 2 Vijk ~ 2S r Svijk + S 2 v 
The first two terms are bounded as before by (|7.57j) and (|7.76|) . For the third we use the 
equation (|7.61| for Vijk, which gives 

(7.77) l^fcl^rVKfcl+ISrWyfcl+inVjfcl+r- 1 ]T l^flV^i.al+r^^^^al 

a+b<2 

We bound all the terms on the right via (|7.57[) and (|7.76l) to obtain 

(7.78) \S^v ijk \<l, i + 4j + 16fc<m-40. 

Now we repeat the I = 2 argument to prove by induction that 

(7.79) \S l r v ijk \<l, i + ij + 16k<m-20l. 

We remark that at each application of the equation P T v = g, one or two S r factors are 
replaced by either an f2 factor or by rr or simply by 1. 

(ii) Again the case m = follows directly from (|7.57j) . Consider now the case I > 1. 
For large r we compute 

rd T {ve" r ) = Sv e lTr + r(d r - ir)ve lTr 
Hence if I = 1 then we need to show that 

(7.80) \(d r +ir)v\ < M" 1 ^ 2 , r > It]' 1 

Repeating the computation at the second and higher levels, the conclusion follows if we 
prove the more general bound 

(7.81) \(d r +iT) l v ijk \ < iTj-V -1 -', r > Irl" 1 , i + 4j + 16fc + 20Z < m 

For I = 1 we rewrite the equation for v in the form 

(7.82) (d r — ir){d r + ir)(rv) = r _1 A w w + rQi r v + rQ sr v + rg 
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By the outgoing radiation condition (|7.58[) we can view this as an ode for (d r +iT){rv) and 
integrate from infinity. By (|7.57[) we can bound the first three RHS terms by |r| _1 r -2 
and integrate to obtain (|T|f) _1 as desired. For the g term a further integration by parts 
is needed. Precisely, we have 

/•OO ■ 

rg{r,uj)e-" r dr = - l -T- l r Q g{r ,oj)e^ Tra + / -T^re^ d r {g{r, oj)e lrT )dr 

r Z Jro 1 

Hence by (|7.T5[) we obtain the (|r|r) _1 bound, and (|7.80p is proved. We further observe 
that we can apply the same argument to Vijk using the equation (|7.61j) to obtain 

(7.83) \{dr+ir)v ijk \ <r" 2 , r > M" 1 

For I = 2 we apply another d r — ir operator in (|7.82[) . This yields 

(d 2 + T 2 )(d r + ir){rv) = r^ 2 A UJ (d r + ir)v + rQi r (d r + ir)v + rQ sr (d r + ir)v 

+ Qi r v + Q sr v + {d r + ir)(rg) 

In view of (|7.57j) and (|7.83j) we can bound all terms on the right by |"r| -1 r -3 except for 
the last one, and obtain (|7.81|) for I = 2 integrating from infinity. For the last term, an 
additional integration by parts is needed, exactly as above. Applying the same argument 
to v^k we also obtain (|7.8ip for / = 2. The bound for higher I is proved in an inductive 
manner. □ 



8. Conclusion 



By Lemma l3.7[ the time Fourier transform u of the solution u to the Cauchy problem 
is represented as 

(8.84) fi(r) =R t (tu + 2P 1 u + u 1 ) := R T (rf + g), f G Z m+1 >\ geZ m ' 2 . 

The bound (|1.12|) will be obtained from pointwise bounds on u(r) and its derivatives with 
respect to r by reverting the Fourier transform. This requires slightly different arguments 
for small and for large r. 

8.1. The case of large r, |r| > 1. If r is large, then a good low frequency approximation 
for u is u f« t -1 /. Subtracting that, we have 

u(t) = T^ 1 / + Ul 

where 

P tUi = (P 1 / + g) + r- x P 2 ! := h + T- 1 g 1 
with /i G H rn ' 2 and G H m ~ 1:3 . Reiterating we obtain the representation 

K-i 

(8.85) u = J2 r-*-Vfc + r- K R T (Tf K + 9k) 

where f — f, fk G H m+1 ~ k>k+1 , and g k G H m+1 ~ k ^ k+2 . Here we choose K large but 
with K <C m. This leads to a decomposition of the high time frequencies P>i(|£)t|)u in 
u of the form 

K-l 

P>i(|A|)« = Ua + ut, u a = x>i(|r|) T ~ k ~ X h, Ub = X>i{k\)T~ K Rr(rf K +g K ). 

k=0 
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For the summands wc have < (r) 3 k therefore applying an inverse Fourier transform 
we obtain 

(8-86) K| < (t)- N (r)- 3 

For the resolvent expression in Ub we use Proposition l7.17l We obtain 

\{rd T )\u b {r)e-)\<\T\- K+l {r)- 1 - 1 . 
Then, inverting again the Fourier transform, wc obtain the following bound: 

(8.87) \u b \<(r}- 1 (t-r) 1 - K 

Together, ([5^6]) and (fST57|) give the high frequency part of ([TTT2"]) . 

To bound the high frequency part of dtu the same argument applies with one exception, 
namely the first term fx>i( T ) arising in the expression for dtu a . However, such a term 
is natural since u is extended by to negative times therefore dtu contains a fSt=o 
component. Factoring that out we obtain as desired 

(8.88) |P >1 (|A|)(5 t u a - fS t=0 )\ < (t)- N (r)-\ 
i.e. the high frequency part of (|1.13[) . 

8.2. The case of small t, |t| <C 1. Here we will still need the arguments from the large 
r case, but there is an additional layer of the proof, given by the resolvent expansion 
around r — provided by Proposition ^. 151 Since / £ ^ m + 1 > 1 and g £ Z m ' 2 , by applying 
Proposition 16.1 51 fa) for / and Proposition 16.1 51 (b) for g wc obtain the representation 

u(t) = (v + TVl + Te(r,T))e lrr + R T ( X> \ T \-i(Tf + g)) + T 2 R T h 

with h as in f|6 . 54[) . Correspondingly we decompose 

P<i(\D t \)u = u a + u b + u c + u d 

where 

Ua(r) = X<l(l r l)-Rr(x>|r|-i(T/ + .9)), = ^ X< 1 ( I A )Rt h 

u c {t) = x<i(|t|)(w + T Vl )e lrr , u d {r) = TX<i(M)e(r, r))e lrr 
For the first term u a we use an expansion as in (|8.85|) to write 

Ua = X<l(M) T " fc_1 X>|r|-i/fc + T- K X<l(\T\)R T (x>\ T \-i(TfK +9k)) 

k=0 

Since fk £ jjm+i-k,k+i an( j g k g jjm+i-k,k+2 ^ ^ e seconc [ term above can be included 
in Ub- We remark that the spatial cutoff X>|r|- 1 plays an essential role, as it allows the 
spatial decay of fx and gx to compensate for the t~ k factor. Hence we redenote 

K-l 

Ua = X<l(M) ^ T ~ k ~ 1 X>\r\-ifk 
fc=0 

Using the pointwise bounds \fk\ < (r)~ 3 ~ fe we invert the Fourier transform to obtain as 
in the large t case 

N / / \ \ N 



< M" 3 



r ) ( 7TT~I I > l^a-P<l(|A|)A= 



r)+tj \ (r) + t 
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For Ub we use Proposition 17.181 to obtain 

|(r9 T )^ b |<r 2 , M^r)- 1 

and 

\{rd T )\u b e^)\<rr-\ |r|>(r)- x 

which lead to 

M < (r)-\t - r)- 2 , \d t u b \ < (r)- x (f - r)~ 3 
For u c we have the pointwise bounds \vq\ < 1, \v\\ < 1 therefore we easily get 

\u c \ + \d t u c \<{r)- 1 {t-r)- N 

Finally we consider Ud, which in this analysis represents the leading term of the con- 
tribution to u of the long range terms in the operator P. 

We recall that the expression e(r, r) has the form 

e(r,r) = (r) _1 ei(r A |r| _1 ) + r(e 2 (r A It) -1 ) - eadrp 1 )), ei,e 2 € S'(logr) 
We represent ei and e 2 in the form 

ei2(r) = ^2c k x>k(r), \c k \ < 1, 

fc>0 

where x>fe are bump functions which equal for r -C 2 m , equal 1 for r ^ 2 m and are 
smooth on the 2 m scale. We denote by e k the contributions of each summand to e and 
consider three cases depending on the size of r: 

i) r « 2*. Then e k has the form 

e k (r,r) = c k T 2 x <2 -" (M) 
therefore its corresponding output in Ud satisfies 

\u k d \ + 2 k \d t u k d \ < 2- 3fe (l + 2 k t)- N 

which leads to (|1.12[) and f|l . 13[) after summation with respect to k with 2 fc > r. 

ii) r ^> 2 k . Then e k has the form 

e fe (r,r) = c fc r(r)~ 1 x< 2 -'= (|r|) 

therefore its corresponding output in Ud satisfies 

\u k d \ + 2 k \d t u k d \ < {r)~ 1 2- 2k {l + 2 k t)- N 

which leads again to (|1.12[) and (| 1 . 1 3|) since now we are summing with respect to k for 
2 k < r. 

iii) r sa 2 k . Here we get terms as in both (i) and (ii), but there is no summation with 
respect to k. 

The proof of Theorem II .41 is concluded. 
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